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OBSERVATIONS OF INTERPLANETARY ENERGETIC CHARGED PARTICLES 
FROM GAMMA-RAY LINE SOLAR FLARES 

M.E. Pesses^, B. Klecker2, G. Gloeckler^, and D. Hovestadt2 

1. University of Maryland, College Park, Maryland 20742 USA 

2. Max-Planck-Institut fur Extraterr. Physik, 8046 Garching, FRG 

ABSTRACT 

Observations of interplanetary energetic ions 
(H+; 5 5 Eh i 20 MeV, He, C, 0, Fe; 1 ^ Ei ^ 20 
MeV/nucleon) and electrons (115 5 Ee ^ 1300 keV) 
from the 7 June, 21 June and 1 July 1980 v-ray 
line solar flares are presented. The observations 
are from the Max-Planck-Institut/University of 
Maryland Ultra Low Energy Wide Angle Telescope 
aboard the ISEE-3 spacecraft. Both June flares 
produced relatively low intensity proton events 
at earth with peak intensities at 10-20 MeV 
-v-5 X 10”2 protons (cm^ sec sr MeV)‘l. Neither 
flare showed evidence of being enriched in either 
^He or Fe at '^'l MeV/nucleon. The 1 July flare 
produced no observable ion or electron enhance- 
ments. 

1. Introduction 


Gamma-ray line emission from the sun was first observed in the great 
flares of 4 August and 7 August 1972 by Chupp et Solar flares which 

produce y.'&y line emission provide a unique opportunity to study ion 
acceleration phenomena in flares. For only in such flares can information 
on the accelerated ion population within the solar atmosphere be obtained 
by direct observations .2 This information confined with interplanetary 
energetic charged particle data can provide the most stringent constraints 
on ion acceleration phenomena. 

The purpose of this brief paper is to present the initial results 
from the Max-Planck-Institut (MP I) /University of Maryland (UM) experiments 
on ISEE-3 on interplanetary energetic charged particles from the 7 June, 

21 June and 1 July y-ray line solar flares (GRLF). 

2. Instrument Description 

The energetic particle data discussed in this paper are from the 
MPI/UM Ultra Low Energy Wide Angle Telescope (ULEWAT), aboard the ISEE-3 
spacecraft. The ULEWAT is a double dE/dx versus E, thin-window, flow- 
through proportional counter/solid state detector composition telescope. 
For a detailed description of ULEWAT see Hovestadt et al_.3 

3. Gamma-Ray Line Flare Identification 

From the launch of the Solar Maximum Mission (SMM) spacecraft, 14 
February 1980, through the period for which we have reduced our ULEWAT 
data, 20 July 1980, the SMM Gamma Ray Experiment (GRE) detected 
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lines from three flares. The photon characteristics of these 6RLF are 
summarized in Table 1. 


T<bl» 1. GanM'Ray L4n« riare CtiaracterUtIct 


Fljr« Date 
(1980) 

Ho 

onset 

X-ray 

Na* 

Location 

■Actiye Classification 
Region Ha X-ray 

Radio OurstsH) 

TT m IV 

Sh* ♦ 2h aJ.y MeV Intens1ty(?)(3) 
7^ inttqral/Aug. 73 

7 June 

0309 

0312 

N14W70 

2495 

SB 

M7 

V 

Y 

Y 

■vl %10‘* 

21 June 

0048 

0119 

N17H91 

2502 

18 

X2 

V 

Y 

Y 

,l«) 

1 Julyf^ 

1618 

1626 

S12W38 

2544 

18 

X2 

Y 

Y 

N 

l1/2 ■v5k10‘* 


(1) V refer* to observed, N not obierve'^. , , , , 

(2) Value* given are for tlie ratio of >eak [integral] intensity of the ■. 'h * y{2.2 MeV) 
line to that observed in the 4 August 1972 uRlf. 

(3) > 0*1 Pyan « al 

(4) 'ha peak of the 2.2 NeV Hr - wv obscured bv rentinuum vrays. When the line emerged 
from coritinuum the ir.tensit ne 1 July i960 GRLF. 

(5) this GPII wet also a white light flare. ‘ 


We have identified the interplanetary energetic particle event de- 
picted in Figure 1 as originating from the 7 June 1980 GRLF. This identi- 
fication is based on:(l) a time delay of 21 i 1 minutes between the hard 
x-ray maximum at the sun and the arrival time of 0.44-1. 3 MeV electrons at 
ISEE-3,(2) the W70° heliographic longitude of the GRLF, and (3) the absence 
of any other large hard x-ray flares within two hours of the observed par- 
ticle event. No solar wind data at earth are available for this period. 

Figure 1 presents 3 
hour averages of several 
ULEWAT counting rates. The 
solid circles represent 
Z > 2 ions with 12 E] ^ 

40 MeV/nuc. The triangle 
represents He with 5 - E^ 
i 20 MeV. The short hori- 
zontal line in the lower 
right hand corner gives 
the upper- limit for He 
5-20 MeV/nuc and Z > 2 
"42-40 MeV/nuc particles. 

We have tentatively 
identified the first of tv/o 
interplanetary energetic 
particle enhancements de- 
picted in Figure 2 as ori- 
ginating from the 21 June 
1980 GRLF. The onset of 
the second enhancement is 
Indicated by the factor of 
two Increase in 115-300 KeV electrons and 10-20 MeV protons intensity at 
^0 LIT on 23 June. The onset for the Z ' 2 {^0.5-12 MeV/nuc) ions in this 
enhancement is not as clear and probably occurs ^^3-6 UT on 23 June. The 
identification of the first event in Figure 2 with the 21 June 1980 GRLF 
is based on; (1) a time delay of 34 i 1 minutes between the hard x-ray mini- 
mum at the sun and the arrival of 0.44-1.3 MeV electrons at ISEE-3, and (2) 
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Figure 2 


the W91° heliographic longitude of 
the GRLF. No solar wind data, at 
earth, are available for this period. 

We are unable to find any in- 
terplanetary energetic particle 
event that is associated with the 
1 July GRLF, in eUher the ULEWAT 
data or data from the UM/MPI instru- 
strument aboard the earth-orbiting 
IMP-8 spacecraft. No solar wind 
data, at earth, is available for 
1 July 1980. However, there are 
solar wind data from IMP-8 for 30 
June 1980. The ^‘'served solar wind 
velocity, on this day was - 300 
kms"l.° This c'^rrosponds to an 
earth-sun flux tube connection at 
-W77®while the GRLF occurred at 
W38“, a difference of 40*^ in connec- 
tion longitude. 


4. Interplanetary Fluxes and Energy S pectra 


An analysis of the GRE '-ray data from the 3 GRLF indicates that pro- 
tons were accelerated up to at least 10-20 MeV^. While it is not possible 
here to make a detailed analyses of the flux, energy spectra and anisotropy 
time histories of the 7 June and 21 June particle events in order to esti- 
mate the injected particle population at the sun, we will present particle 
fluxes and spectral indices at the time of maximum particle intensities. 
This information is given in Table 2. The values here are derived from 3 
hour average counting rates. The >'s are calculated from a least squares 
fit to the equation J = where *- = Enerqy/nucleon , and J is the differ- 

tial intensity. For the 1 July 1980 GRLF only upper limits of the particle 
flux can be given. 




Table 

Peak Particle FI 

u« and Speitra Indices 


riar** 

n.ite 

P^oton^ (fir 

."C STr Ili'Vl' 

■ ric; 

trrns (rT" se. strl‘1 

rroton Llectron 

t > 

S-lO 

10-.’l» MrV 

115-100 kek 

300-140 KeV 0.44.l.TMrV 

? Jui'r 

1 5.10*' 

6.B.10*’ 

rs.ui’ 

31.10’ 3.M0' 

: 1 1.5 







ll 

3.1.10 ‘ 

1 (»>io 

6. .’.lO*^ 

I s.m‘ z.r.io' 

IB ?,6 

I Ju1> 

.1.6.10*’ 

1 :.io*’ 

r.5 

7 6.10'’ -1 4.10'’ 
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5. E lemental an d Is otopic Composition _ 

Because of the low Z - 1 ion flux in both 7 June and 21 June 1980 
GRLF associated particle event we are able to determine composition ratios 
at equal energy per nucleon in only one ULEWAT energy range. 

The relative composition ratios are presented in Table 3. 
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For the 7 June event the 
data are summed from 7 
June 1200 to 8 June 2359 
For the 21 June period 
the data are accumulated 
for only the first event 
in Figure 1 from 21 June 
0000 UT to 22 June 1200 
UT. For comparison we 
have also presented com- 
position ratios from the 
August 1972 GRLFs, for 
^He rich flares, normal 
composition flares, the 
photosphere, and the 
corona . 

6. Dis cussion and Con - 
?H7sions 

We conclude from 
results presented in 
Table 3 that neither the 7 June 1930, 21 June 1980 nor the August 1972 
GRLF are classifiable as ^He rich or as Fe rich. This is consistent with 
the plasma instability preheating model of the ^He and Fe rich events 
discussed by Fisk^. We ,‘lso note that the GRLFs in Table 3 compared to 
normal composition flares i re depleted in Fe relative to 0 by a factor of 
2-3. 

Although the intensity of ,-ray line emission at 2.2 MeV from the 
1980 GRLFs is the same order of magnitude as that from the August 1972 
GRLF, the peak vlO MeV proton intensities in the two spts of flares differ 
by many orders of magnitude. Furthermore we have observed several non 
GRLF associated particle events with ULEWAT in which the neak proton in- 
tensity exceeds that of the GRLF associated events by 2 orders of magni- 
tude. If the proton flux difference between events is not due to propa- 
gation effects it could ifidicate that there is little i^ any relationship 
between the ion population that give rise to the %-ray line emission and 
the ions population that is injected into interpl anetary rpace. 
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TIME AND ENERGY DEPENDENCE OF HEAVY ION ABUNDANCES 
IN SOLAR FLARE ENERGETIC PARTICLE EVENTS 

G.M. Mason^, H. Welss^, G. Gloeckler^, and D. Hovestadt^ 

iDepartment of Physics and Astronomy, University of Maryland 
College Park, MD 20742 USA 

2Max-Planck-lnst1tut fur Physik und Astrophysik 
8046 Garching b. Munchen, W. Germany 

ABSTRACT 

We present examples of the time and energy dependence 
of the abundances and spectra of the major heavy Ions 
He, C, 0 and Fe during solar flare events, taken from 
a survey using the UHD/MPI ULET telescope on IMP-8 
during 1973-1977. In some cases, we find time variations 
in the 0/He, 0/C and Fc/0 ratios which appear to be in- 
consistent with models based solely on rigidity dependent 
propagation in the interplanetary medium. 

1. Introduction 

In recent years, the availability of relatively low threshold (^-1- 
few MeV/nuc) heavy-ion detectors have made possible the first exploratory 
investigations of heavy-ion abundance variations during individual solar 
flare events (Armstrong ^ al. 1976; O'Gallagher et 1976; Scholar et 
a^. 1978; vonRosenvinge and Reames 1979). A coirmon, although not univer- 
sal feature identified in these early studies has been a decrease in the 
Fe/0 ratio (taken at the same kinetic energy/nucleon) following the flare 
onset, which can be modeled in terms of rigidity dependent propagation in 
the interplanetary medium when taking account of the partially stripped 
Fe vs. the (nearly) fully stripped 0. However, other abundance variation 
features have been identified in isolated flare events which do not appear 
to be as readily understandable, e.g., 0/He variations (vonRosenvinge and 
Reames 1979; Weiss and Mason 1980), and spatial variations (Armstrong et 
al. 1976). 

To form an overview of these abundance variations, we have carried 
out a survey of 13 major solar particle events during the recent solar 
minimum (1973-77) using the ULET telescope on IMP-8 (described by Klecker 
et 1977). Several of the flare periods such as July 3-10, 1974, Nov. 
S^IO, 1974, and Sept. 10-18, 1977, exhibited complex variations of the 
.nonitored abundance ratios whose interpretation is not obvious, and which 
might be caused by some coirtoination of source characteristics along with 
propagation and shock event effects in the interplanetary medium. The 
most commonly (though not always) observed trend in the periods is a ten- 
dency for ratios of (heavier ion/lighter ion) to decrease as time passes; 
the rate and magnitude of this decrease depends on the ratio examined. 
Below we give examples of two relatively simple events which show this 
behavior. Even in these cases, though, the details of the abundance var- 
iations are difficult to understand by considering only the differing 
rigidities of the particles. 

2. Observations for the 20-23 September 1974 Flare 

The 20-23 Sept. 1974 particle event was apparently associated with a 
IB flare in McMath region 13225 (N08W63) at 2253 on Sept. 19. The upper 
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Fig. 1 . Count rates, spectral indices and abundance ratios during 
19-28 Septefnbtr 1974 

left panel of Fig. 1 shows ULET 3-hour average count rates for Z ^ 1 and 
Z i 3 ions of ^0.5 MeV/nucleon for this event as well as a smaller one 
starting 23 Sept, (not discussed here). A shock wave, marked S in the 
fig., passed IMP-8 at 1244 on Sept. 21 (M. Pesses, private comnunication) . 
The lower left panels of the figure show 3-hr. average spectral indices > 
for the monitored fluxes over the range of 0.6-2. 5 MeV/nuc, assuming power 
laws in kinetic energy/nucleon: dJ/dE = (const)E''*. The spectral indices 

show clear velocity dispersion effects, which are not unexpected for a 
flare of this type, which occurred in the so-called well connected region. 
Notice however, that major changes in the Fe spectral index > occur more 
quickly than for the lighter species, and it is larger at later times than 
the indices for the lighter species, leading to a decrease in the Fe/0 
ratio at higher energies which may be easily ' een in the plotted abundance 
ratios. For He, C and 0, the changes in appear to stop after the shock 
passage. 

The right panels of Fio. 1 show 3-nr. average abundance ratios of 
(heavier species/lighter species) during this period for each of the 3 
energy windows monitored. The figure format for abundance ratios arranges 
kinetic energy/nuc windows vertically, and ratio types horizontally result- 
ing in 3 repetitions of the time scale. The Fe/0 ratio in all energy win- 
dows is seen to decrease rapidly (1/e fall -off time Tq 12 hrs) by a fac- 
tor of -vB after the flare onset, consistent with rigidity dependent propa- 
gation of partially stripped Fe ions as discussed by O'Gallagher et al . 
(1976) and Scholer ^ (1978) for this event. After this raplT^faTl 

off, the time variatTon becomes less rapid. In the case of the 0/C ratio. 
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there is also a decrease, at least up to the time of the shock passage, 
with fall -off time To 100 hrjrs. Tinally, the 0/He ratio is also seen 
to systematically decrease in time witn Tq ''' 40 hrs, and in this case the 
rate of decrease does not appear to be changed by the shock passage. 

Observations for 22-26 Noventer 1977 

A second example of time dependences of the type discussed above is 
seen in a November 1977 particle event associated with a 2B flare at 0945 
on 22 Nov. in the "well-connected" region at N24W40. A sudden comence- 
ment was observed at Earth at 1227 on 25 Nov. Fig. 2 shows the count 
rates, spectral indices and abundance ratios during this period. Velocity 
dispersion effects are clearly evident in the spectral index data from the 
onset time until about 1200 on 23 Nov., after which the spectra continue 
to steepen, but at a reduced rate. As in the precedin’, case, the Fe/0 
ratio quickly drops after the onset (Tq -v 15 hrs) and then becomes more 
nearly constant. It is interesting to note in this case that the change 
in temporal behavior of the Fe/0 ratio is not associated with the shock 
passage, while in the case of the Sept. 1974 event such an association 
might be inferreo. Up to the ti\?e of the shock passage, the 0/C and 0/He 
ratios decrease in a similar manner to the Sept. 1974 events. There ap- 
pears to be a definite change in 0/C associated with the shock passage. 

We emphasize that the time variations in the abundance ratios in this 
event are not ordered by the individual spectral indices. That is, the 
time decreases in the ratios do not appear to depend upon the heavier ion 
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having a steeper (or flatter) spectral slope than the lighter Ion, since 
It can be seen that Fe/0 decreases with time with >Fe > ''^o* 0/^* 

decreases with time with "Tq < 

4. Discussion 

The clear velocity dispersion effects seen In these events suggest 
an Interpretation based on propagation effects In the Interplanetary medi- 
um, where at the same velocity partially stripped Fe has a hlaher rigidity 
than nearly stripped 0 {O'Gallagher et al. 1976; Scholar et 1978). In 
order to examine the new data presenfedlTere In this light, we use the 
most recent observations for solar flare Ion charge states for HeV/nuc- 
leon particles, from Hovestadt ^ al. (1900, 1981) who find 
mean charge states of He, C, 0, anT^e of 1.9, 5. 8, 7. 2 and 13.5, respective- 
ly. With these charge states, the ratio of particle rigidities at equal 
velocities are P(0/He)«l .06; P(0/C)-1.07; P(Fe/0)»l .87. Thus, the signif- 
icantly higher rigidity of Fe vs. 0 might well lead to more rapid propaga- 
tion for Fe in a model where the scattering mean free path Increases with 
rigidity, as described in earlier works. However, the 0/He and 0/C varia- 
tions are puzzling considering the small rigidity differences of the 
species, as is the fact that the 0/He variations are more similar to the 
Fe/0 than the 0/C variations, while consideration of the rigidity ratios 
would lead to the opposite expectation. We conclude that rigidity differ- 
ences alone cannot explain the observations presented here. 

We speculate that other factors, such as abundance differences at the 
flare site on the sun, or features of the acceleration and release mech- 
anisms may play an important role in the variations discussed here, as 
well as the mere complex cases in flares cited In si. We have also found 
that in general the individual flare spectra are very often not fitted 
well by the kinetic energy/nucleon power laws used here: not only are 

the spectral indices for the species often different from each other -- 
leading to abundance variations at differing energies -- but also the 
spectra tend to "roll over" towards lower energies. We intend to examine 
other spectral forms (e.g. distribution functions in particle energy/ 
charge) in order to search for a simpler representation of the data, and 
to help separate the numerous different effects operating in these parti- 
cle events. 
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MEASUREMENTS OF IONIZATION STATES FOR 
LOW ENERGY HEAVY IONS IN SOLAR FLARE PARTICLE EVENTS 

L.S. Ma Sungl, G. Gloeckler^ and D. Hovestadt^ 

1. University of Maryland, College Park, MO 20742 USA 

2. Max-Planck-Institut fur Extraterr. Physik, 0-8046 Garching FRG 

ABSTRACT 

We examine the distributions of ionization states for flare associa- 
ted heavy pjrticles (Z > 2) between 76 and 1000 keV/charge. Tht particle 
data have been collected by the University of Maryl and/Max-Pl ancL- 
Institut Electrostatic Energy-Charge Analyzer (EECA) Experiment on IMP-8. 
The EECA sensor measures particles in the same energy/charge interval. 

By examining the particle energy pulse-height distribution, the ioniza- 
tion states of Z ' 2 ions can be deduced if the elemental abundances are 
known. We determine the elemental abundances between 1.0-4. 6 MeV/nucleon 
from particle data collected by the Ultra Low Energy Telescope (ULET) on 
the same experiment. In this energy ranae, a power law fit to the parti- 
cle spectra would give energy dependent abundance ratios. However, if 
the particle distribution function f is plotted as a function of the 
quantity n, where n = vRi^, v and R are the particle speed and rigidity 
respectively, and n is an adjustable parameter, all elements are found 
to share the same distribution f % e’n/no with a com.Tion no ^nd nH) for 
a given particle event. In this repres*‘ntation , the abundance ratios 
become independent of the spectral shape and can therefore be applied 
to the lower energy range cf the EECA data. For the ionization-state 
distribution, we assume tfiat ions are in charge equilibrium at a given 
coronal temperature(^) . By comparing our data in several energy-per- 
charge windows with calculated distributions which are produced by 
folding the instrumental profiles with an incident particle population 
of given elemental abundances and a charne-state distribution, we can 
determine the charge distribution characteristic of the particle event 
and hence deduce the temperature of the source materials and the possible 
acceleration process for these ions. 

1. Gloeckler, G. , Weiss, H., Hovestadt, D. , Ipavich, F.M.^ Klecker, B., 
Fis^ . L.^*., Scholer, M. , Fan, C.Y. and O'Gallagher, J.J. , 1981, 17th 
Inte.-nat ional Co^smj_c_ Ra^_Cqnf. (Paris, France). SH 3.1-7. 

2. Jordan, C. 1^9, Mon. Not. R. Astr. Soc., 142, 501-521. 
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OBSERVATIONS OF THE IONIZATION STATES OF ENERGETIC PARTICLES 
ACCELERATED IN SOLAR FLARES 

6. Gloeckler^, H. Weissl, D. Hovestadt^, F.M. Ipavich^, B. Klecker^, 
L.A. Fisk3, M. ScholerS C.Y. Fan^, and J.J. O’Gallagher^ 

1. University of Maryland, College Park, MD 20742 USA 

2. Max-Planck-Institut fur Extraterr. Physik, D-8046 Garching, FRG 

3. University of New Hampshire, Durham, NH 03824 USA 

4. University of Arizona, Tucson, AZ 85721 USA 

5. University of Chicago, Chicago, IL 60637 USA 

ABSTRACT 


We report results of a survey of ten solar flare particle 
events in which we measured trie ionization states and spec- 
tra ot 0.3 to 2.4 MeV/nuc He, C, 0 and Fe. We find He* to 
be present in all events (He'^’/He’*"'’ 'v 0.1 to 0.25), and the 
mean ionization states of C, 0 and Fe to be high. The dis- 
tribution functions of He'*’, He’*"*' and heavier elerrients are 
well represented by simple exponentials of the particle 
speed times its rigidity to a power n, where n is between 0 
and 1, and equal e-folding values. Our results are consis- 
tent with a model whereby ions are accelerated in the corona 
by multi-dimensional shocks out of a population taken from 
both hot (Tg /?, 2* 10° K) and cold(Tg lO^ K) coronal regions. 


1. Introduct ion 

Knowledge of the ionization states of particles accelerated in solar 
flares is important r>jt only because it provides information about the 
characteristics of the regions where acceleration takes place but also 
because it is essential to know the charge to mass ratio of these parti- 
cles in order to determine the nature of the acceleration and propagation 
mechanisms. Earlier measurements of charge states of -t.1 MeV/charge solar 
flare accelerated ions (e.g. Gloeckler et Ma Sung et ^.3) indi- 

cated high ionization states, consistent with a coronal material at a 
temperature of 1-2-106 K. The recent discovery by Hovestadt of a 

surprisingly la'^ge relative abundance of He'*' at higher energies, on the 
other hand, requires a much cooler source material. 

In this brief report we outline results of a survey of 10 solar 
flare particle events observed over a one year period starting in Sept. 
1978, using the ULEZEQ sensor of the Max-Planck-Institut/University of 
Maryland experiment on ISEE-3 positioned vO.Ol AU upstream of the earth. 

We find that energetic he* is present in all of the flare events examined, 
being the third most abundant element after H and He’*"*'. The mean charge 
states of C, 0 and Fe, on the other hand, are consistently high, the aver- 
age values over all 10 flares being 5.8, 7.1 and 13.5 respectively. 


In strumentation 

In the ULEZEQ sensor the chemical elements are first identified 
using pulse-height analyzed signals from the proportional counter (dE/dx 
element) and the solid state detector (E element) of the dE/dx vs E tele- 
scope placed at the exit of an electrostatic deflection system wFich dis- 
perses ions according to their energy per charge. The charge states of a 
particular element are then determined from the amount of deflection 
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Ionization state histoprams for He, carbon, oxygen and iron for 
the June 6-9, 1979 solar flare particle event. 

measured by a position sensitive detector, and the residual energy E. 

This analysis is performed for each incoming ion. The instrument and 
analysis technique are more fully described by Hovestadt et al 

3. Obse rva tions 

In order to illustrate the results provided by our survey we discuss 
in some detail particle observations for the June 6-9, 1979 event which 
are typical of all the flare events examined. Fig. 1 shows the charge 
state histograms for helium, carbon, oxygen and iron in the indicated 
energy per nucleon intervals. The resolution of the instrument allows us 
to separate the individual charge states of He although not of the heavier 
elements. The He'*’ abundance is about 10« that of He'*’^ . The ionization 
states of the heavier elements are high; the mean charge states of C, 0 
and Fe are 5.7, 7 and 12.5 respectively. Our preliminary analysis does 
not allow us to estimate the distribution of charge states of each heavy 
element, although it is evident, especially in the case of iron, that 
several ionization states are present (e.g. 10 to 15 for Fe). However, 
we find no evidence in any of the heavy ion charge histograms examined 
for low ionization states (e.g. O'*’, Fe^^) above a few percent level. 
Finally, we point out that our present determination of mean ionization 
states is subject to systematic errors of .^.05 of the charge state. 

The energy spectra of He'*"*', He'*, 0, C and Fe are shown for the June 
1979 event in Fig, 2(a). In the restricted energy range of 0.3 to 2.4 
MeV/nucleon the spectra may be fitted to power-laws, j = AE"'». We note, 
however, that the spectral indices are different (e.g. 1.7 for He'*"*’ 

and 2.2 for Fe) and that the power-law fit is especially poor for He'*'*. 

In order to find a spectral representation in which the energy dependence 
of the ratios is removed we show in Fig. 2(b) the dependence of the dis- 
tribution functions, f = j/v^, on the quantity vR*'* where v and R are the 
speed and rigidity of the ions respectively, and n is a parameter which 
is adjusted to minimize the energy dependence of the ion abundance ratios. 
With n = 1/3 for the June event we find excellent fits to all the data, 
where the distribution function for each species i has the form 


n 
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fi » Aiexp(-vR^/Eo) with a conmon Eq “ 0.26(MeV/nuc)2/3, jp computing 
the rigidities of 0, C and Fe we used mean ionization states derived from 
the charge histograms of Fig. 1. 

A sunt.iary of our results for the 10 flare particle events is given 
in Table 1 which lists the relevant flare characteristics, spectral infor- 
mation, HeVHe'*"*' ratios and mean ionization states for C, 0 and Fe. In 
the power-law representation of the differential intensities, j, the spec- 
tra of He'*’ are consistently steeper than those of with a flare aver- 
aged ratio of the spectral indices ''^He'^’/^He'*"*’ 1.17. In a representa- 

tion where the distribution functions f « expC-vRi^/Eo) . however, one can 
find values of n for which Eq is conmon to He'*’ and He'*^'*’. Typical values 
for both n and Eq are ^.4. Although the He''’/He'*’'*’ ratio is energy depen- 
dent when computed at equal energy per nucleon, this ratio is constant in 
the representation where the species have a common distribution function. 
Assuming the vR*’ dependence of the distribution function is valid over a 
broader energy range (i.e. 0.1 ^ E 10 MeV/nuc) we can derive density 
ratios for He^/He'*’^ by integrating f over velocity space. These ratios, 
listed in column 10, are seen to range from M3. 08 to 0.29 with a flare 
average of 0.15. Finally, the mean ionization states of C, 0 and Fe, de- 
rived from charge histograms such as shown in Fig. 2, are found to be con- 
sistently high, in general agreement with previous measurements at lower 
energies (e.g. Gloeckler ^ ^.1*2, Ma Sung et al .3) . 

4. Discussion and Conclusion 

The fact that energetic He'*’ is seen at surprisingly large relative 
abundances in every solar flare particle event examined indicates its com- 
mon presence in the source material from which solar particles are accel- 
erated. One plausible explanation of this result and of our finding that 

Table 1 Spectral Cbaracterittics and lonltation States of 0.3 to 2.4 f4eV/nuc 
Me, C, 0 and In Solar Flare Particle Ivents 


Flare 1 

Characteristics 

^Pdctrel 

Characterist ics 



loniratlon States 


(b) (b> 

Me Me Ne*^ 

n 

*0 

Me’/Me** 

Chiroe SUles 

Date 

laportance 

Location 

(M 
■^1 'A 

ric) 

n 

rid) 

C 

c(f) 

0(f) 

Fe**’ 

1978 













2S-27 Sept 

38 

wso 

1 S4 1 

05 

0.50 

0.55 

0.06 

0.22 

0 11 

5.7 

7.2 

13 

10-12 HOv 

ll“’ 

H38 

1 a: 1 . 

.15 

0.53 

0.33 

0.06 

0.17 

0.08 

6.0 

7.2 

14 

1979 













l»-22 Fib 

In'-’ 

U.’7 

2 36 1 

06 

0.38 

0.33 

0.06 

0.18 

0.10 

5.7 

7.3 

15 

2»-29 Mar 

in'“’ 

IS6 

2 25 1 

16 

0.4S 

0 37 

0.13 

0 55 

0.29 

5.6 

6.9 

-- 

3-6 Apr 

ib‘*' 

U14 

1 . 78 

23 

0 47 

0.47 

0.06 

0.17 

0.09 

5.8 

7.2 

13 

23- 2S Apr 

7 


2 S*- I 

27 

0 71 

0 48 

0 09 

0.50 

0.21 

6.F; 

7.2 

-- 

28-30 May 

Ih'*’ 

all 

2 42 1 

05 

-0.05 

0 IS 

0.12 

0.09 

0.10 

5.9 

7.2 

-- 

6-9 June 

26 

Mb 

1 66 1 

21 

0.35 

0.42 

0.08 

0.25 

0.14<*’’ 5.7 

7.0 

12.5 

10 July 

Id**’ 

£sr 

1 72 I 

25 

0 42 

0.44 

0 10 

0.33 

0.16 

5.7 

6.9 

13 

1&-N l«t 

3 

1107 

LL? L 

JO 

0.60 

q_B2 

0.13 

0.50 

0.23 

5.5 

LI 

— 

FLAM AflRAfit 



1.92 1 

17 

0 44 

0.44 

0.09 

0.30 

0 15 

5.8 

7.1 

13.5 


(a) HiUlpla flare Injections 

(b) anerv dependant, avcraMd over 0 3 to 2 4 Neil/nucleon 

i c) ratio at edual n • vd" (v • velocity, R • rigidity) 

dj dinsity ratio attualng f • Aeap(-n/(o) Molds froe [ t 0 ■ to t { 10 NeV/nuc 

t) poitibly entrqjr depen^t, ratio averaged over 0 4 to 2.4 MaV/nucleon, values subject to tystflaatlc 
error of 0.05«Q». 

(f) one slyaa error Is ‘V.IS charge units 
(tj one t1^ error Is 'vl.O charge units 
(h) density ratios He*^;Me* 0 C Fe are SS 71 0.6 0.16 
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JUSE 6-9, 1979 
ULEZEO SENSOR on ISEE*3 
Moi * PIvKk * Inctiiuli/ 
Uhivtrsity of Mvytond Eipcrincnt 

>5 



(MeV/Nuc) 


10 2C 30 40 50 
(\^OCi»y) «(Rpdity)° 


vR 


0 55 




_2. (a) Energy spectra and (b) 
bution functions of He"^, 

C and Fe in the June 6-9 


distrl- 
He-*-, 0, 

19 79 event 


the 


the distribution functions 
f depend on particle rig- 
idity, is that ions are 
accelerated by a shock gen- 
erated as a result of the 
flare eruption. This shock 
accelerates material not 
only from the hot coronal 
regions (which contain He'*"*’, 
0+6-8, Fe"*^!®"!^, etc.) but 
also from the much colder 
regions (Tg % 10^ K) where 
He^, 0"*^ etc. would be pre- 
sent. Cold regions have in- 
deed been observed (Bohlin, 
private communication), and 
f « exp(-vR^/Eo) is consis- 
tent with acceleration by 
multi-dimensional shocks 
(Lee, private communication). 
The reason we do not observe 
low charge states for the 
heavy elements is that in the 
limited energy range of our 
measurements the abundance 
of singly ionized heavy par- 
ticles, which have a high 
rigidity, would be signifi- 
cantly reduced relative to 
particles. For example, if we 


lower rigidity, highly ionized heavy 
assume for the June 1979 event that f is the same for O'* as it is for 
0’*7.4, apd that the density ratio of "cold" to "hot" material is 0.14 
(see column 10 of Table 1), we find that only v7 O'* ions would fall within 
our energy range. This represents less than 0.5% of the total oxygen 
population and would be beyond the detection limit of our instrument. 

This clearly demonstrates the importance of measuring ionization states 
over a broader energy range with better charge resolution and sensitivity. 
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EVOLUTION OF THE ENERGETIC PARTICLE COMPOSITION DURING 
SOLAR FLARE EVENTS AS OBSERVED BY VOYAGER 1 AND 2 

D.C. Hamilton and G. Gloeckler 

Department of Physics and Astronomy, University of Maryland 
College Park, MD 207^2 USA 

ABSTRACT 

We present measurements of the composition and energy spectra of the 
energetic ions from two solar flare particle events and discuss their 
evolution with tine during each event. The data are from the Low Energy 
Particle Telescope (LFPT) on Voyager 2, which can identify the major ion 
species (Z = 1-26) over an energy range of '0.5-50 MeV/nucleon. 

A prompt particle event resulting from a 2B flare at N23W40 on 22 
November 1C77 was observed by Voyager 2 at 1.6 AU frcn the sun. The 
event averaged ene *gy spectra (covering 6'r days) of H, He, C, N, 0 and 
Fe, for differential flux expressed as a power law in kinetic energy per 
nucleon, typically have spectral indices in the range 1. 5-2.0 from 0.6 
to \2 MeV/nuc and steepen at higher energies. Abundance ratios evaluated 
at equal energy/nucleon (0.60-0.95 MeV/nuc and 2. 1-3.2 !1eV/nuc) display 
the following characteristics; the 0/He, 0/C, and Fe/0 ratios all have 
maximum values at event onset and decrease fairly smoothly throughout 
the event, the amount of decrease depending on species and energy. These 
results are qualitatively similar to those reached from observations at 
1 AU in a similar energy rjnce'^'. The behavior of the He/H ratio, 
hcv/ever, is more complicated, exhibiting a minimum value ('-.014 for 0.60- 

0. 95 MeV/nuc) at onset, increasing to a maximum value ( .067) at time of 
maximum flux, and then gradually decreasing during the decay phase to a 
value of "1.029 late in the event. 

A large solan particle event was observed by Voyager 2 in September 
1979 at ^5.6 AU from the sun. This particle event does not appear to be 
associated with an optical flare and most likely was produced by a flare 
on 14 Sept, in McMath region 12698 which was beyond the east limb of the 
sun at that time (at E107 ). An v? x-ray event peaked at 0802 on 14 
Sept, and had associated Type II and IV radio bursts. The particle 
event at 5.6 AU was characterized by hard energy spectra and long dura- 
tion. Event averaged spectra for H, He, C, N, 0, and Fe (covering 17 
days) typically have spectral indices in the ranne 1.3 - 1.7 over an 
energy range of 0.6 to 20 f‘eV/nuc, steepening only above 20 MeV/nuc. 

In general, the onset portion of the event is complicated, the rise to 
maxirum intensity being dominated by spatial structure rather than IP 
propagation. Voyager 2 was apparently not well connected to the flare 
site. The low energy particles (‘^a few MeV/nuc) show evidence of 
enhancement due to passage of at least one IP shock (possibly associated 
with a CIR) during the event onset. After time of maximum intensity, 
however, durino the relatively sm.ooth decay phase of the event, almost 
all the abundance ratios of heavier to lighter species (He/H, 0/He, 

0/C, Fe '0) at energies from 0.6 - 1.8 MeV/nue decrease slowly from their 
values at maximui'. intensity to values up to a factor of 2 smaller late 
in the event. 

1. Mason et al.. 1981, 17 ICPC Paris, paper 3.1-1. 
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ON THE ANT I CORRELATION BETWEEN THE ^He/^He RATIO 
AND PROTON INTENSITY IN 3He RICH FLARES 

M.E. Pesses 

Department of Physics and Astronomy, University of Maryland 
College Park, Maryland 20742 USA 

ABSTRACT 

A model is presented which can explain the observed 
anticorrelation between the size of the ^He/^He 
abundance ratio at ^10 MeV/nuc and the peak intensity 
of '40 MeV protons in ^He rich events. If ^He is 
preaccelerated with respect to ^He prior to the main 
ion acceleration phase, and if the main ion acceler- 
ation phase has ar effective injection cutoff which 
varies from flare to flare, the observed anticorre- 
lation follows naturally. 

1. Introduction 


Solar flare energetic particle events in which the abundance of ^He 
relative to ^He is enhanced by several orders of magnitude over coronal 
and photospheric abundances have been the subject of much discussion. 1 *2 
A particularly interesting feature of the ^He rich solar flares is the 
observed occurrence of the largest relative enhancements of 3ne to AHe 
in the events in which jjhe peak proton intensity at ^-10 MeV is the lowest 
protons (cm^ sec sr MeV)~l. On the other hand, no 
3he rich flares have been observed in events in which the peak proton 
intensity is large ' lO^-A protons (cm2 sec sr MeV)"l. 

The relationship between the intensity ratio of ^He to ^He at equal 
energy per nucleon rpHe/^He), and the peak event intensity of '>-10 MeV 
protons Ip, has been studied by Hempe et ^.,3 and Ramaty ejt al.2 The 
find that over the range of 2xl0"2 rl^He/AHe) 10^ and, 10^ Ip ^ 

1q 2 protons (cm2 sec sr MeV)"^ the observations can be fit by a 

power law, r(JHe/AHe) = constant x Ip"'^ where > '>- S* 

The purpose of this paper is to show that this anticorrelation is a 
natural consequence of producing abundance enha^'cement by a two step 
acceleration process. 

2. The Model 

The 3ne rich events have been interoreted as being due to the nuclear 
spallation of '10 MeV/nuc “^He into and due to the selective 

heating of 3He by plasma waves prior to the main acceleration phase® 

The nuclear reaction model is apparently unable to explain either the 
observed nonenhancement 2 h and 3 h or the observed enhancement of Fe and 
other heavy ions relative to 0 in 3he rich flares. 

For this reason I am assuming in this paper that the 3 hc enhancements 
are due to a selective preacceleration of 3 hc prior to the main ion accel- 
eration phase. The main ion acceleration phase accelerates ions to the 
energies at which they are observed in interplanetary space. 

A schematic representation of the relative number density of ^H, ^He 
and 3He prior to the main ion acceleration phase is shown in Figure 1. 
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The solid lines are the number densities 
oT 1 h, 3h 6 and *^He as a function of some 
particle variable, like velocity, or 
rigidity or energy/charg®, before pre- 
acceleration occurs. The wavy line is 
the number density of ^He after the pre- 
acceleration occurs. 

The relationship between the number 
density at the main acceleration site 
prior to the main acceleration phase n^ 
and the observed number density of parti- 
cles at 1 AU Oq can be expressed by 
using the acceleration-transport parti- 
cle propagator P, which gives the proba- 
bility of a particle with some initial 
value of a variable x between x-j and 
xi+dxi being observed at a value of the 
variable between Xq and XQ+dxQ. The re- 
lationship is given in equation (1) for 
the case where x = energy/nuc = T. 

no(To)«To = h‘dnj(T,)/dT,P(To.Ti)'iTidTo 

, ( 1 ) 

where nolTjJdTg [ni(Ti)dT.j] is the obser- 
ved [initial J number density of particles 
with observed (initial] energy/nuc between 
Tq and Tp+dTo [Ti and Ti+df^]. The rela- 
tionship between the observed integral intensity from Tq = T^ to Tq = T 2 , 
1^2 ^>^d no(Tp)dTp is given by Equation (2), 

I 12 = (2To/mp)'^^(To)dTo (2) 

where mn is the proton mass and the particle observed velocity »2To/mp 
is << than the soeed of light. 

As pointed out by Fisk^, orders of magnitude enhancements in ^He re- 
lative to ^He can only be produced if there is some cutoff on particle 
injection into the main ion acceleration phase. If particles of all ini- 
tial velocitv, rigidity , energy per charge etc., haye an equal probability 
of being accelerated to Tq there will be no large -^He enhancements as the 
initial and final values of r(^He/^He) will be basically the same. There 
should be some differences in the abundance ratio as most acceleration 
processes are, in theory, dependent in some fashion on particle mass to 
charge ratio, which is 3/2 for and 2 for 

When the injection cutoff is applied to protons as well as helium 
the observed ant 1 correlat ion between ri^He/'^He) and Ip can be shown to 
follow naturally. . ^ 

Assume for alaebraic simplicity that P(To,T,) = ( '-^UTg (Tj/To)^ for 
T,- Tq and zero -‘or T, • Tc and T, ■■ Tq, where T^ is the injection cutoff 
energy 'nuc, and i is a constant ^ 0. Also assume for algebraic simplicity 
that dni(T.j)/dTi = K(Ti./T*)"', for T, > T^-, where T* is a constant, and 
K and '• do not depeno on T ^ or T* but can depend on species. Then Equa- 
tion il 'becomes 



s' V ^ 

\ * 

\ 


N . 
F" 

Figure 1 
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"o(To)<JTo ■ K(T,/T*)-''(a*l)T;'(T,/T„)“dT,dTo 

As Tq 10 MeV/nuc, is very likely >>T^., which for y 
to 

no(To)dT„ 1 K(a-U)(i-l-a)''(Tc/T*)-'(VT^)‘'*'dTo 


(3) 

reduces (3) 

(4) 


Combining (4) with (2) gives, 

1^2 [2^,5^ (T2-Tj)(a-‘i){>-l-a)]'^(T‘p-T^'")K(a+l)T*> (5) 


then r(^He/4He) » Ii 2 i^or ^He) ♦ Ij 2 (^or ^He) which gives 

r(3He/^He) ^ K3K^^(>4-l-ci)(>3-l-a)'^(T^/TM’'4‘^3 (6) 

And also Ip = Ij 2 (^or ^H), which gives 

( 7 ) 

In (6) and (7) the subscripts 1, 3, and 4 on K and > refer to the values 
of K and > for iH, ^He and ^He respectively after preacceleration and 
prior to the main ion acceleration phase. In Equation (6) > y^- 


3. Results and Discussion 

Equations (6) and (7) show that as T^- increases in value r(^He/^He) 
increases in value while Ip decreases in size (if > 1 + a), which is 
consistent with the observed relationship between and Ip dis- 

cussed in Section 1. This result can be visualized from Figure 1. As 
the acceleration cutoff varies in value the number of ^He and 4He ions 
which are injected into the accelerap’on process also varies. Because ^He 
has been preaccel erated relative to ^He, it has harder spectra. Hence 
as the cutoff velocity, rigidity, energy/charge etc. increases the number 
of ^He relative to 4He that are injected into the main acceleration pro- 
cess increases, while the number of injected into the main acceleration 
pnase decreases. Since n^ increases as n^ increases the observed anti- 
correlation between r(3He/4He) and Ip follows naturally provided T^. varies 
from flare to flare. For example, for one value of Tr in Figure 1 
r(“He/4He 10"^ and for a larger value of Tj. rl^He/^He) 1, with the 
size of Ip decreasing as the size of ’'(“He/4He) increases. 
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ANALYTICAL DESCRIPTION OF CHARGED PARTICLE TRANSPORT ALONG 
ARBITRARY GUIDING-FIELD CONFIGURATIONS 

J.A. Earl 

Department of Physics and Astronomy, University of Maryland 
College Park, MD 20742 U.S.A. 

ABSTRACT 

A new desrripticn w*” focused transport has been developed 
which is valid for arbitrary spatial dependences of both the 
scattering mean free path and the focusing length. In partic- 
ular, it describes the "supercoherent transition" from the 
coherent modes that occur when focus' ng is strong compared 
to scattering to the diffusive mode that occurs when focusing 
is weak. Because it specifies the pitch angle distribution 
in terms of relatively simple functions, the new formulation 
facilitates not only the detailed interoretation of solar 
particle events, but also the qualitative understanding of 
focused transport. In the steady state, it predicts that 
the angular distribution consists of a highly collimated 
"strahl" which is superimposed on an isotropic "halo". This 
picture provides new insight into observed configurations of 
steady-state particle fluxes. 

I. Introduction . Although particles are clearly scattered by random 
interplanetary maonetic fields, the mean free path V during many solar par 
tide events is comparable not only to the distance to the sun, but also 
to the focusing length L which describes spatial variations of the guiding 
field B, and whicn is defined by 

(1/L) = -(1/5)(3B/ -2). (1) 

wnere z is distance parallel to B. To obtain a better description of the 
non-diffusive effects that occur under these circumstances, we have anal- 
yzed particle transport under the assumption that both the scattering mean 
free path > and the characteristic length L for guiding field variations 
are constant in space and time. (Earl 1976a. 1976b-. Bieber 1977) These 
assumptions do not accurately reflect conditions in interplanetary space, 
but they drastically simplify the mathematics of the transport problem. 
This simplified formulation has explained many features of solar particle 
events (Ma Sung and Earl 1978, Bieber ^ 1979, 1980), but some discre- 

pancies have appeared in the detailed correspondence between theory and 
experiment, "hese discrepancies clearly indicate the need ^or an improved 
formulation. 

In response to this need. Ng and Wong (1979, 1981) and Gombosi and 
Owens (1980) calculated solar particle fluxes with L dependent upon z. 

The method introduced by these workers is a valuable contribution, but 
its purely numerical ^orm encumbers the physical understanding of its im- 
plications. Consequently, this paper attempts to provide the insight 
inherent in an analytic approach by formulating the problem of focused 
transport in terms of partial differential equations that 'pply to 
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arbitrary spatial variations of \ and L. The present analysis follows 
lines similar to those already taken by Kunstmann (l‘J79), but it differs 
from his treatment, because it considers two components instead of one, 
it examines the coupling between these modes, and it includes the effects 
of spatial variations of \ and L. 


2. The New Transport Equations. Particle transport is described, 
in terms of* the distribution function f, by the Boltzmann equation, which 
can be written in two alternative forms. 


+ yV — = 4 e^ — e"^ — 
at ^ az ? 3u au 


3(f/B) , ,, 3(f/B) . 1 3 r, (2) 

3t 3i: 


where the odd function G is given by 

G{z,u) = -(V/B)('^B/3z) (l-t'2)(dr/Mi/}), (3) 

0 

and where t is time, V is particle velocity, u is the pitch angle cosine, 
and ;{u,Z( is the Fokker-Pl anck coefficient of pitch angle scattering. 


In the calculation of f, it makes no difference which form of the 
Boltzmann equation is solved, for they are mathematically equivalent. 
However, certain integrals over the two forms give completely different 
results. More specifically, the integral over m of the second form gives 
the familiar flux equation, which states that the temporal rate at which 
the density changes is equal to the negative divergence of the flux. If 
the first form is multiplied by exp{-G} and integrated, the result ex- 
presses a similar but less familiar relationship between a weighted den- 
sity and a weighted flux. Both of these independent relationships must 
be satisfied by any valid solution of the Boltzmann equation. 


A first-order approximation to the distribution function is: 
f{p,z,t, = F,-z,f’-^BH.{z,t}e^, (4) 


In the steady state with * and L constant, which represents the zeroth- 
order situation, the coefficients Fr and H.- are independent of space, 
time, and each oiher. Equation (4) describes a slightly perturbed situ- 
ation in which both coefficients are weakly dependent upon z and t, ard 
in which v and L are weakly dependent on z. Under these circumstances, 
the coefficients Fc and Hr, cannot be chosen arbitrarily, for the distri- 
bution function f must be consistent with the integrals over the Boltz- 
mann equation that were discussed above. This requirement leads to the 
relationships. 
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These arc transport equations, whose solutions give a useful description 
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of particle transport, and in which D is the coefficient of spatial diffu- 
sion. Their parameters are given by 

K = S e^’ du, and = (V/2K) ue^ du. (7) 

Equations (5) and (6) give the familiar phenomenon of diffusion in 
the weak focusing limit, (\/L)<<l, where \V/3L)=( D/L) , but they give 
essentially coherent transport when focusing is strong, (x/L)>l, and 
V,i^V. This evolution in structure unfolds continuously as X and L vary 
smoothly with z. This means that solutions of the first order transport 
equations automatically take into account the changes in the qualitative 
nature of particle transport that occur in physically realistic models of 
interplanetary conditions. At the Munich conference, for example, I 
described these changes in terms of a "supercoherent transition" at which 
the nearly coherent propagation of flare particles in regions of strong 
focusing near the sun gives way to ordinary diffusion in the outer solar 
system (Earl 1976b). The new formulation makes it possible to deal with 
this transition on a quantitative basis. 

4. Steady State Solutions . It is not possible to find analytical 
solutions of the general time-dependent first-order equations. At this 
stage, numerical metiiods are appropriate, but they have not yet been im- 
plemented. Instead, a simple, new, and general solution of the steady 
state problem is described below. 

In the steady state, the temporal derivatives in equation (6) vanish, 
and its solution takes the form, 

H. = (j/V„K), (8) 

where a is a constant. With this result, equation (5) takes a form, 

B.' D 9B ,g) 

■ ■ (V.K? B ’ 

which can readily be integrated to give expressions for the isotropic 
component . 

The figure shows the results of such a calculation for a monopolar 
guiding field configuration, B a r"2, in which v decreases as ' i r"2/3, 
and into which particles are steadily injected at r = 0. Because the 
ratio ( '/!) = l- decreases with the coordinate r, which replaces z, this 
configuration embodies a smooth transition from coherent to diffusive 
transport, "he densities, Fr , BKH^ and Fq, are all plotted as ratios to 
the density Fq predicted under the’same circumstances by the diffusion 
equation. Here, the total density Fq is the sum of the isotropic den- 
sity Fr and BKH^ , which is the density embodied in the anisotropic com- 
ponent. The abscissa is the ratio of r to the radius r^ at which 
(v/L) = 1. 

Far from the origin, where ('/L)--! and transport is essentially 
diffusive, B<H^ is small compared to Fr and their sum Fg is identical to 
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the diffusive result. At 
the top, a polar plot of the 
angular distribution p = 1 
exemplifies the relatively 
weak anisotropies that occur 
in this region. Near the 
origin, where (r/ro) < 0.1, 
transport is nearly coherent, 
the anisotropic component 
dominates, and its density 
is much larger than that 
predicted by diffusion 
theory. The strong aniso- 
tropies which appear in 
this region, and which are 
illustrated by the angular 
distribution for p = 10, 
correspond to a highly 
collimated beam of particles 

^ '0| 01 ' / 1 10 10 streaming away from the 

source. 

In the transition region, 0.1 < (r/ro) density components 

are comparable, and the angular distributions, as illustrated by the plot 
for p = 4, consist of a collimated "strahl" superimposed on an isotropic 
"halo". These are exactly the features observed aboard HELIOS in the 
angular distribution of solar wind electrons by Rosenbauer et al^. (1977), 
who interpreted these features in terms of unscattered and backscattered 
populations of electrons. The present analysis puts this interpreta- 
tion on a solid analytic foundation, and opens up the possibility of 
deriving quantitative information from observacions of steady-state angu- 
lar distributions. 
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A SUPRA-THERMAL ENERGETIC PARTICLE DETECTOR (STEP) 

FOR COMPOSITION MEASUREMENTS IN THE RANGE 
''^20 keV/nucleon to I MeV/nutleon 

G.M. Mason and G. Gloeckler 

Department of Physics and Astronomy, University of Maryland 
College Park, Maryland 20742 USA 

ABSTRACT 

A novel detector system is dest/ibed, employing a time- 
of-flight versus residual energy technique which allows 
measurement of particle composition (H-Fe), ener^ spec- 
tra and anisotropies in an energy range unaccessible 
with previously flown sensors. Applications of this 
method to measurements of the solar wind ion composition 
are also discussed. 

1. Introduction 

In studies of particle acceleration and transport in astrophysical 
plasmas, phemonena occurring within our solar system may offer the best 
hopes of understanding many of the responsible physical mechanisms, since 
these sources are available for close and detailed observations with numer- 
ous information channels. Unlike the cosmic ray accelerators, however, 
these local sources accelerate particles to modest energies — usually 
less than a few MeV/nucleon. Composition measurements of these low energy 
particles help identify the particle sources; while details of the compo- 
sition, energy spectra and anisotropies make it possible to distinguish 
amongst possible models for acceleration and transport. Thus, low energy 
composition measurements are essential in the study of many particle accel- 
eration sites including, e.g., the solar atmosphere, planetary magneto- 
spheres, interplanetary shock waves, and co-rotating interaction regions. 

Previous instruments designed for detailed composition measurements 
at low energies have employed the dE/dx vs. E technique, wherein the dE/dx 
detector element was kept as thin as possible in order to achieve a mini- 
mum energy threshold. The limits of present technology appear to have 
been reached with thresholds for, e.g.. He of 'v-lSO keV/nucleon (Krimigis 
et a^. 1977; Hovestadt et 1979), and with good resolution of ion 
species above a few hundred keV/nucleon. In this paper we describe a new 
technique using measurements of time-of-fl ight versus energy (Gloeckler 
1977; Gloeckler and Hsieh 1979), with which considerably lower energy 
particles may be studied. 

2. Instrumental Technique 

Fig. 1 is a cross section of the Supra-Thermal Energetic Particle 
Detector (STEP), which is one of a pair currently under construction for 
flight on the Solar Polar Mission. The thin foils and Si detector in the 
telescope are rectangular (1.5 x 4.0 cm), with the figure showing a cut 
through the 1.5 cm dimension. The inner foil and solid state detector 
(SSD) are separated by 10 cm, and define a f ield-of-view (17° x 44°) with 
a geometrical factor of 0.4 cm^-sr. 

Ions entering the system pass through the foils and stop in the low- 
noise SSD. The signal Emeas the SSD is related to the ion energy 
deposit E by Emeas = “E, where a is the "defect" in the SSD signal re- 
sulting from non-ionizing collisions of the ion with the Si lattice. The 
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value of a is determined in labora- 
tory calibrations (Ipavich et al . 
1978). 

In addition to the energy de- 
posit, the ion velocity is measured 
by making use of secondary electrons 
which are emitted with a few eV when 
the ion passes through the surface 
of the inner foil, and the front of 
the SSD. The secondary electrons 
are accelerated to 1 kV, focused, 
and deflected to hit the surfaces of 
microchannel plates (MCP) positioned 
to intercept electrons from the foil 
("Start" MCP) or the SSD ("Stop" 
MCP). The electrostatic accelera- 
tion and focusing are done using 
techniques developed in vacuum-tube 
technology (e.g. Harman 1953), and 
avoid the use of accelerating grids. 
The MCPs produce large, fast signals 
suitable for determining the ion's 
time-of-flight, t, using appropriate 
fast electronics. 

Having measured the particle 
energy deposit, Emeas/“» ^nd time- 
of-flight (TOP) over a known dis- 


( 1 ) 


The atomic number Z is not determined; however, for all major species over 
the range H-Fe, a mass measurement is sufficient to identify the ion uni- 
quely. The incident energy of the ion is obtained after taking account 
of the ion's energy loss in the foils. With this information, plus the 
known pointing direction of the spacecraft, the composition, energy spec- 
tra and anisotropies of the incident ions may be determined. 


3. Energy Range and Mass Resolution 

The energy range of the STEP telescope for several species is shown 
in Fig. 2. The lower energy threshold is determined by the SSD discrimi- 
nator threshold (^-30 keV) and foil thickness. For the STEP telescope, 
relatively thick foils (■vlBOO ^ each) are required to limit the sensitiv- 
ity of the system to solar photons. For applications where sunlight can- 
not strike the aperture, a very thin foil can be used, yielding consider- 
ably lower energy thresholds (see §4). The upper limits of. the energy 
ranges shown in Fig. 2 are set by the threshold value of the time-of- 
flight electronics (^3 nsec) for H and He, and by saturation in the SSD 
amplifiers for heavier species. 

It can be seen from equation (1) that the mass resolution of the sys- 
tem is affected by dispersions in Emeas* f STEP, the major 

factor below MeV/nucleon is the dispersion in Emeas* resulting from 
non-ionizing collisions of the ion with the Si lattice. In Univ. of Mary- 
land SSDs similar to those used in STEP, the quantity (FWHM of Etneas)/Emeas 
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Increases from a few per- 
cent at 100 koV/nucleon 
to 20-30% near the instru- 
ment threshold (Ipavich 
and Galvin 1980). Taking 
this effect into account, 
along with dispersions in 
T and d, the instrument 
is calculated to have re- 
gimes of isotope, element 
and element group resolu- 
tion as shown in Fig. 2. 

Fig. 3 shows the 
measured mass resolution 
of the STEP prototype 
telescope, obtained dur- 
ing its first calibration 
with heavy ions. For the 
upper panel, which is 
typical of energies near 
At energies near 200- 


300 keV/nucleon (lower panel) op^/m covers the range from 0.05 ( ^He) to 
0.08 (C). These values were about twice the expected result, due to a 
crosstalk between the MCPs and the SSD which ^doubled the FWHM of the 
energy signal in the SSD over the values measured with the same detector 
when removed from the telescope. It is expected that removal of the cross- 
talk problem in the flight telescopes will result in the resolution shown 
in Fig. 2. 

The large energy range covered by the STEP sensor, and the mass reso- 
lution already achieved in the prototype telescope, show that many hither- 
to inaccessible regions of the low energy charged particle spectra may be 
studied for the first time with 
this instrument. 
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A . Oth er App l i cat ions 

The technique used in STEP 
may be extended down to solar wino 
energies, and may include ion 
charge state determination, by 
adding features shown in Fig. 4, 
which is a cross section of the 
Solar Wind Ion Composition Spec- 
trometer currently under con- 
struction for flight on the Solar 
Polar Mission (Gloeckler ejt 
1977). In SWICS, solar wind ions 
of energy E, mass m, and charge 
state q enter a collimator and 
pass through a deflection system 

whose voltage is stepped through 
a cycle to select particles of 
different E/q. The particles are 
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then accelerated by a 
voltage Va(- 30 kV) be- 
fore passing into the 
TOP system. Sunlight 
is blocked by the col- 
limator, making it pos- 
sible to use an extreme- 
ly thin (”^5 ug/cm) mesh- 
supported carbon foil 
in the TOP system. 

After passing through 
the foil, the particles 
. . - have an energy IqVaI + 

V E', where E' is the 
• -* * •*■ — * - original energy less 

the small energy loss 
of the ion in the foil. 
As in the case of STEP, 

Pig . 4. Cross section of Solar Wind Ion Composi- ion's flight time 
tion Spectrometer (SWICS). ^ distance d is 

measured, as well as 

the output (Emeas) of the SSD. The mass of the ion is then given by 
equation (1). The ion charge state q is obtained from 

q = (En)eas'a)(!VAl ^ E'/q). (2) 

The incident energy E is found using the value of q determined in (2), 
along with the E/q setting of the deflection system. With this informa- 
tion, and the known pointing direction of the spacecraft, the ion's charge 
state, energy, mass, and arrival direction are determined. 
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